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Abstract. We present here the results of Arecibo timing of PSR B1516+02B, a 7.95-ms pulsar in a binary system with a 
~ 0.17M,;, companion and an orbital period of 6.85 days located in the globular cluster M5. The eccentricity of the orbit 
(e = 0.14) has allowed a measurement of the rate of advance of periastron: a> = (0.0136 ± 0.0007)°yr . It is very likely 
that the periastron advance is due to the effects of general relativity; the total mass of the binary system is (2.14 ± O.16)Af . 
The small measured mass function implies, in a statistical sense, that a very large fraction of this total mass is contained in 
the pulsar: M p = (1.94+j} [g)M Q (1-rj); there is a 5% probability that the mass of this object is below 1.59 M®. With the 
possible exception of PSR J1748— 2021B, this is the largest neutron star mass measured to date. When combined with similar 
measurements made previously for Terzan 5 I and J, we can exclude, in a statistical sense, the "soft" equations of state for 
dense neutron matter, implying that matter at the center of a neutron star is highly incompressible. There is also some evidence 
for a bimodal distribution of MSP masses, the reasons for that are not clear. 
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TIMING OF M5B 

Over the past 20 years, more than 130 pulsars have been 
discovered in globular clusters (GCsU Among the first 
discoveries were PSR B 1 5 1 6+02A and PSR B 1 5 1 6+02B 
1 1]. Both of these millisecond pulsars (MSPs) are located 
in the GC NGC 5904, also known as M5; for this reason 
we will refer to them as M5A and M5B. The latter is 
a 7.95-ms pulsar in a binary system with a ~ 0.17 — Mq 
companion and an orbital period of 6.86 days. At the time 
of its discovery, this was the MSP with the most eccentric 
orbit known (e = 0. 14). In the Galaxy, 80% of MSPs are 
found to be in binary systems, and with a single excep- 
tion (PSR J 1903+0327, see David Champion's contribu- 
tion to these Proceedings) they are in low-eccentricity or- 
bits with white dwarf (WD) companions. In GCs, grav- 
itational interactions with neighboring stars can greatly 
increase the eccentricity of these binary systems |2j]; the 
orbital eccentricity of M5B is ~ 10 4 — 10 5 times larger 
than that of Galactic MSP- WD systems with similar or- 
bital periods. When Anderson et al. (1997) published the 
timing solutions of M5A and B, they used the eccentric- 
ity of M5B to detect its periastron advance, but the large 
relative uncertainty of the measurement did not allow any 
astrophysically useful constraints on the total mass of the 
binary. 

In yj, Freire et al. report the results of recent (2001 
to 2007) L-band observations of these two pulsars. The 
first 2001 observations were part of an Arecibo search 



1 See Scott Ransom's review, in these Proceedings. For an updated list, 

see |http : //www2 . naic ■ edu/~pf relre/GCpsr .html| 



for pulsars in GCs, which found a total of 13 new MSPs 
J5I. Three of these were found in M5, subsequent obser- 
vations of this GC were made chiefly with the aim of tim- 
ing the new discoveries. However, both M5A and M5B 
are in the same radio beam as the new pulsars. They are 
clearly detectable in the L-band data, allowing for tim- 
ing of much better (M5A) or comparable (M5B) quality 
to that obtained at 430 MHz. Including those previous 
data, this provides a much longer total timing baseline 
(18 years) and much improved timing parameters. 



Observations, data reduction and timing 

The L-band observations started in 2001 June. Until 
2003, we used the "old L-Wide" receiver in the Gregorian 
Dome (T sys = 40K at 1400 MHz). Since 2003 February, 
we have been using the current "L-Wide" (T sys — 25K at 
1400 MHz). The Wide-band Arecibo Pulsar Processors 
(WAPPs, la]) made a 3 -level digitization of the voltages 
of a 100 MHz-wide band for both (linear) polarizations, 
correlating them for a total of 256 lags. These were then 
integrated for a total of 64 jxs and the results of both 
polarizations added in quadrature and written to disk. 
At first, only one WAPP was available. In this case we 
centered the observing band at 1425 MHz. After 2003, 
three more WAPPs became available, and we started 
observing simultaneously at 1170, 1420 and 1520 MHz, 
thanks to the wide frequency coverage of the new L-wide 
receiver. 

The lags were then Fourier transformed to generate 
power spectra. These were then dedispersed at the known 
DM of these pulsars and folded modulo their spin pe- 



riods using the PRESTO pulsar software package^. We 
then cross-correlated the resulting pulse profiles with the 
average pulse profile in the Fourier domain [7] to ob- 
tain topocentric times of arrival (TOAs). These were then 
analyzed with TEMPC0, together with the TOAs derived 
from the old 430-MHz observations made from 1989 to 
1994. We used the DE 405 Solar System ephemeris H to 
model the motion of the Arecibo 305-m Radio Telescope 
relative to the Solar System Barycenter. 

The resulting timing parameters are presented in J3]. 
Astrophysically, the most important new parameter is the 
rate of advance of periastron for M5B: ri) = (0.0136 ± 
0.0007 )°yr- 1 . 



Binary, pulsar and companion masses 

As argued in Q], the CO is solely due to the effects of 
general relativity. In this case, we can estimate the total 
mass of a binary system: 
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where T Q = GM /c 3 = 4.925490947^is. Using the ri) 
above, we obtain M = (2.14 ± O.16)M . For the nomi- 
nal ri) and a median inclination of 60°, the mass of the 
companion is O.166M and the mass of the pulsar is 
1.98M . This is well above all the neutron star masses 
that have been precisely measured to date. 

We calculated a 2-D probability distribution function 
(pdf) for the mass of the pulsar and the mass of the com- 
panion, assuming that the pdf for ri) is a Gaussian with 
the half-width equal to the 1-C7 uncertainty of the ob- 
served cb and an a priori constant probability for cos i. 
The two-dimensional pdf is then projected in both di- 
mensions, resulting in 1-D pdfs for the mass of the 
pulsar and the mass of the companion. These are dis- 
played graphically in Fig. Q] The pulsar definitely has 
a mass smaller than 2.35M , and the companion has 
a mass larger than 0.13M Q , the median and 1-c lim- 
its for the pulsar and companion mass are 1.94+g JZjMq 
and O.164+J]^ 2 M respectively. There is a 99%, 95% 
and 90% probability that the pulsar is more massive than 
1.19, 1.59 and 1.69 M . There is a 1.3% probability that 
i is low enough to make the neutron star mass fall within 
the range of NS masses observed in double neutron star 
(DNS) systems: from 1.20 Mn measured for the com- 
panion of PSR J1756-2251 B to 1.44 M measured 
for PSR B1913+16 |[Toh. 



http : / /www . cv . nrao . edu/~s ransom/ presto 
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STATISTICAL EVALUATION OF MASS 
MEASUREMENTS 

A list of mass estimates for millisecond pulsars is pre- 
sented in Table Q] Most of these were derived for MSPs 
in GCs. The reason is, of course, than in GCs the MSP- 
WD orbits can be perturbed, resulting in much higher 
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FIGURE 1. Constraints on the masses of M5B and its companion. The hatched region is excluded by knowledge of the mass 
function and by sinz < 1. The diagonal dashed lines correspond to a total system mass that causes a general-relativistic 0) equal or 
within 1-ct of the measured value. The five solid curves indicate constant inclinations. The gray bars indicate the range of precisely 
measured neutron star masses (from ~ 1.20M Q to 1.44M ft ). We also display the probability density function for the mass of the 
pulsar {top) and the mass of the companion (right), and mark the respective medians with vertical (horizontal) lines. 




eccentricities than possible in the Galaxy. These eccen- 
tricities allow the measurement of post-Keplerian (PK) 
effects like co. Eventually other PK effects will be mea- 
surable, like the Einstein delay (7). However, GC pulsars 
are generally rather faint, therefore such measurements 
will require some time. 

For that reason, and as in the case of M5B, most 
of the estimates in Table Q] are based on measurements 
of cb only. These are "incomplete" measurements, in 
the sense that one more PK parameter is necessary to 
have an unambiguous determination of the pulsar mass. 
Nevertheless, unambiguous upper limits for the mass of 
the pulsar and lower limits for the mass of the companion 
can always be obtained in these cases. Furthermore, in 
systems where the mass function is small (as for M5B) 
there is a much greater probability of most of the mass of 
the binary belonging to the pulsar itself (see above). 

From Table Q] we can see that the assumption that 
the CO is relativistic yields for M5B the largest neu- 
tron star mass presently known, with the possible ex- 
ception of PSR J 1748 -202 IB (NGC 6440B, see Paulo 
Freire's contribution on the new pulsars in NGC 6440 
and NGC 6441 in these Proceedings). 



Statistical evidence for high neutron star 
masses 

One of the interesting features of Table [T] is that as 
the total binary mass increases, the mass function does 
not increase, with a single exception: PSR J05 14— 4002A 
(see Paulo Freire's contribution on NGC 1851A in these 
Proceedings). This system has a highly eccentric orbit 
that is thought to have resulted from an exchange inter- 
action [18]. If the increase in total mass for the other sys- 
tems was due solely to an increase in the mass of the 
companion, then there should be a more general trend to 
higher mass functions, i.e., PSR J0514-4002A should 
be the rule among the massive binaries, not the excep- 
tion. 

It could happen that the separation between the "light" 
and "massive" binaries (shown by the horizontal line in 
Table[TJ is due merely to a significant classical contribu- 
tion to CO in the latter group. This could explain all the 
anomalously high masses produced by the co estimates. 
However, it is unlikely that the extra classical contribu- 
tions have exactly the values required to make the mass 
estimates nearly identical to each other. Furthermore, a 



possible cause of an extra contribution to cb is tidal ef- 
fects due to an extended companion. If the presence of 
extended companions were to explain the large measured 
values of cb, then it should also lead on average to tidal 
circularization of the orbits of the "massive" binaries, 
particularly for those with shorter orbital periods (even 
if that is not guaranteed in individual cases, like M5B, 
which could have been perturbed recently). In reality, the 
"massive" binaries are not less eccentric than the "light" 
binaries. 



Statistical combination of mass pdfs 

The mass pdfs for the 7 MSPs in Table Q] with cb 
measurements were calculated as described above for 
M5B, and displayed graphically in Fig. [2] The mass 
medians for Ter 5 I, J and NGC 6440B are 1 .87, 1 .76 
and 2.73 lfT7h respectively. All these pdfs have long, low- 
probability tails towards the low masses, corresponding 
to improbably low orbital inclinations. 

Combining the pdfs for the masses of Ter 5 I and J, 
lfl5ll reached the conclusion that at least one of these is 
more massive than 1.48 and 1.68M with 99% and 95% 
confidence levels. Combining those pdfs with the mass 
pdf for M5B, the probability that all of these pulsars are 
less massive than 1.72 and 1.79 Mq is only 1 and 5% 
respectively. These limits introduce some of the most 
stringent constraints to date on the bulk behavior of cold, 
super-dense matter. They exclude many of the "soft" 
equations of state that have been proposed to model that 
behavior IU9I1 . 

From NGC 6440B alone, we can derive a 99% proba- 
bility for that pulsar having a mass larger than 2.01 Mq, 
which, if confirmed, would be an even tighter constraint 
to the EOS. If all cos are relativisitc, the probability of 
all these massive binaries having "normal" (1.2 < M p < 
1 .44.1/ . ) masses is 9 x 1CT 9 . 

FORMATION 

The mass pdfs in Fig. [2] suggest that the distribution of 
MSP masses is bi-modal, with NGC 6440B as a super- 
massive outlier. PSR J0514-4002A and the MSPs in the 
"light" (M < 2Mq) binaries have masses smaller than 
1.5Mq, i.e., they are not significantly more massive than 
mildly recycled neutron stars, despite having spin fre- 
quencies of hundreds of Hz. In particular, the case of 
M28 C shows that MSPs can be recycled by accreting 
< 0.15M Q from their companions. Other MSPs are sig- 
nificantly more massive, it is not clear why they are so. 

It could happen that they were born that way. A bi- 
modal (or tri-modal) distribution like that of Fig. [2] is 



exactly what is predicted by hydrodynamical core col- 
lapse simulations [ 20] : stars below ~ 18M Q are expected 
to form ~ 1.20— 1.35 M,., NSs. Stars with masses be- 
tween 18 -20M© form 1.8M NSs. Above 2OM , stars 
experience partial fall-back of material that can signifi- 
cantly increase the mass of the remnant, making it either 
a super-massive NS or, if its mass is above the maximum 
stable neutron star mass, a black hole. 

This possibility raises the question of why such mas- 
sive NSs have never been found in the 9 known DNS 
systems (which have a total of 18 NSs). Given the nar- 
row range of progenitor masses (18 — 20M Q ), it is possi- 
ble that massive NSs are relatively rare. However, 50% of 
the NSs in eccentric binary MSPs in Table[T]are massive. 

It is possible that all these NSs started instead with 
similar masses. In the case of MSPs, the accretion 
episode is much longer than for the recycled pulsars in 
DNS systems, with a potentially (but not necessarily) 
larger mass transfer. This is a natural explanation for why 
we only see massive NSs as MSPs but not in DNS sys- 
tems. If this was the case, we should then expect that 
the more massive MSPs, having accreted more mass and 
angular momentum, should spin faster than the less mas- 
sive MSPs. Table Q] shows that the opposite is true: the 
more massive MSPs spin more slowly (v < 125 Hz) than 
the less massive MSPs (v > 200 Hz). More statistics are 
needed to verify the significance of this relation; but if 
it holds, then there might be two MSP recycling mech- 
anisms, one of them transmitting more angular momen- 
tum and the other more mass. If, on the other hand, neu- 
tron stars start with different masses, and the more mas- 
sive NSs have a higher moment of inertia, then transfer- 
ring the same amount of angular momentum to a massive 
NS will cause a smaller increase in the spin frequency. 
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FIGURE 2. Probability distribution functions (pdfs) for seven of the MSPs in Table Q] The mass pdfs of the MSPs in the least 
massive binaries (those with At < 2M Q ) are represented by the dashed curves. Despite the limitation of being calculated solely from 
0), the pdfs capture well the peak in NS mass that is known to occur at 1.2 — 1.4 Mq, suggesting that the peak at 1.8 — 2.0Mq is a 
also a real feature. 
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